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The innnite-(7 Anderson model coupled to a Jahn- Teller phonon is studied using the slave boson 
method on the basis of the large degeneracy expansion (1/N) scheme. The model Hamiltonian acts 
on the orbital degrees of freedom. The main focus is on the interplay between strong local electron 
correlation and weak Jahn- Teller electron-phonon interaction. The Kondo temperature is found 
to decrease by Jahn- Teller interaction. The influence of the Jahn- Teller interaction on dynamical 
correlation functions is very significant in sharp contrast with the case of the Holstein-type phonon 
which couples to charge degrees of freedom. 
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The orbital degrees of freedom of electrons in a solid 
play a prominent role in modern material sciencei In 
particular, it turns out that electric current can be con- 
trolled by manipulating the orbital states of electrons. 
Spurred by this possibility, a new research area called 
orbitronics is emerging.— The orbital physics is also rel- 
evant to the quantum transport through orbitally active 
atoms (or defects) in nanostructuresi 3 ^ 

The dynamics of the orbital degrees of freedom is de- 
termined by the Coulomb repulsion between orbitally ac- 
tive electrons and the coupling to the vibrations of ions 
surrounding the electrons. A well-known example where 
the orbital physics manifests itself is the transition-metal 
oxidesi Electrons of these materials occupy d orbitals, 
and the Coulomb repulsion between these electrons is 
very strong. At the same time the electrons in d or- 
bitals are coupled to the vibrations of surrounding oxygen 
atoms via the so-called Jahn- Teller (JT) coupling. 5 The 
JT coupling stems from the anisotropy of d orbital wave 
functions. The orbital pseudospin operator to which JT 
phonon couples is defined by 
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where a is a real spin index and a, b are orbital indices. 
For the case of e g orbitals of 3d electrons, the orbital 
indices a,b designate two orbitals d x 2_ y 2 and d 3z 2_ r 2, 
and the matrices r Q are the 2x2 Pauli matrices. 

The progress of theoretical understanding of orbital 
physics has been slow in spite of its prime impor- 
tance. The reason is that the physical phenomena of 
orbital physics are usually associated with strong elec- 
tron correlation and simultaneously with strong coupling 
to phonons of various types. This makes theoretical anal- 
ysis exceedingly difficult. Despite these difficulties, the 
phase diagram of manganese oxides could be qualitatively 
understood by the mean-field approximation study on 
the generalized Hubbard model, 6 but the understanding 
of dynamical properties in low energy strong coupling 
regime is beyond the reach of such mean field approxima- 
tion. In passing we mention that for weakly interacting 



electron systems there exists a very successful Migdal- 
Eliashberg theoryii 

To understand orbital physics in the presence of strong 
correlation we need the approximation schemes which 
preserve the essential features of strong correlations. One 
of such approximation schemes is the dynamical mean- 
field theory (DMFT)&2 In DMFT a lattice Hamiltonian 
is mapped (or approximated) to a certain quantum im- 
purity Hamiltonian. In this brief report we study an An- 
derson impurity Hamiltonian interacting with local JT 
phonons (AJT model) from the perspective of the or- 
bital physics within the DMFT scheme or the transport 
phenomena in nanostructures. 

To simplify our problem further we will assume that 
real spins are polarized, namely, the material under con- 
sideration is in a ferromagnetic state. This simplification 
is not unrealistic, since a large portion of phase diagrams 
of transition-metal oxides and multiferroics are ferromag- 
netic. A typical example is the low temperature phase of 
Lao.7Sro.3Mn03 manganese compound with T c ~ 370 K 
which shows the colossal magnetoresistanceiiS An addi- 
tional advantage of studying the ferromagnetic state is 
that we can focus on the orbital dynamics itself not be- 
ing mixed with spin dynamics. 

The physical properties of the AJT model become 
clearer when they are compared with those of Anderson- 
Holstein (AH) model. 11 ! 12 : 13 ! 14 : 15 ! 16 The dynamical de- 
grees of freedom of the AH model is the real spin. The 
orbital degrees of freedom is not present in AH model. 
Most important difference between two models is that 
the Holstein phonon couples to the charge degrees of 
freedom of impurity electron while JT phonon couples 
to orbital degrees of freedom. It is well known that the 
nonperturbative Kondo singlet ground state of the An- 
derson model (in local moment regime) is due to the 
quantum fluctuations in spin (or orbital pseudospin in 
the context of AJT model) channel. In this regard one 
can expect that the influence of Holstein phonon on the 
Kondo ground state would be small unless the electron- 
phonon coupling is extremely large, and this expectation 
has been confirmed i^i^i^i^i^ On the other hand the in- 
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fluence of JT phonon is expected to be substantial since 
it directly couples to the channel responsible for Kondo 
(orbital pseudospin) singlet ground state. The results of 
our study confirm this expectation. 

The AJT Hamiltonian consists of the following three 
parts: 



(2) 



7i c i and 7i p h are the Hamiltonians for the electron and 
phonon parts, respectively. 7Yjt is the Hamiltonian de- 
scribing the interaction between the impurity electron 
and the JT phonon, 

Wei = - ^) 4a<*a + £ / £ tif* + U fW\h 
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Hj T = 9oQT z , T z = ~(/J/ T - flf L ). 
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The indices a =f , J, denotes the orbital degrees of freedom 
as |= d. 
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and |= c?3 Z 2_ r 2, and the real spins are as- 
sumed to be perfectly polarized. f a is the local impurity 
electron operator, and c ka is the conduction electron op- 
erator. Vk is the hybridization matrix clement, /i c is the 
chemical potential for the conduction electrons, e/ is the 
energy of the impurity level, and U is the local Coulomb 
repulsion at the impurity site. T z is the z component 
of orbital pseudospin. Q is the local Jahn- Teller phonon 
coordinate. P is the conjugate momentum of Q satis- 
fying [Q, P] — ih, and M is ion mass. The JT phonon 
is assumed to be an Einstein phonon with frequency f2. 
iVi a t is the number of lattice sites for the conduction elec- 
trons, go is the electron-phonon coupling constant which 
is assumed to be weak. 

^From Eq. (jSJ) one can notice that the JT phonon plays 
a role of the fluctuating magnetic field acting on orbital 
pseudospin. The coupling to the magnetic field evidently 
hinders the formation of the Kondo singlet ground state, 
thus the Kondo temperature at which the singlet begins 
to form is anticipated to decrease. Our results will show 
that this anticipation is indeed correct. The electron- 
phonon interaction of the AH model is 7Yho1 = 9oQp, 
where p = /J/ a is the impurity charge density, and 
the index a is the real spin index. 

The essential features of the Kondo physics of the An- 
derson model are kept intact in the limit of infinite local 
Coulomb repulsion UmL In this limit the doubly occu- 
pied impurity state is not allowed in the physical Hilbcrt 
space. The elimination of doubly occupied state provides 
a considerable formal simplification in the framework of 
the slave boson method>i£*i2*22i Furthermore it possesses 



the great advantage that it can describe the nonpertur- 
bative Kondo singlet ground state at the mean-field level. 
We study the model Hamiltonian Hajt in the limit of in- 
finite Coulomb repulsion, employing Coleman's slave bo- 
son method in the scheme of 1/N expansioni&SI (N = 2 
for the original 7Yajt)- In this limit the impurity electron 
operator f a can be expressed as 

fl = 4b, (6) 

with a constraint b jf b + J2 a s a s a — 1- b is the slave boson 
operator, and s a is a constrained fermion operator. The 
above constraint is implemented by a complex Lagrange 
multiplier A = zOq + A sa . Next we take the Read-Newns 
gauge in the scheme of 1/N expansion in the form s a = 
z a e te ,b = (^/N/2)re w . The conduction electron and JT 
phonon can be integrated out exactly, and then we are 
led to the following (imaginary time) effective action: 
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where Cf = e/ + A sa , q = 1/N, and 9 = 9 + Q . E is the 
self-energy due to the interaction with conduction elec- 
trons. The imaginary part of So is given by — zAoSgn(e), 
where Ao describes the hybrization amplitude which is 
usually much larger than the Kondo energy scale, while 
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is the bare phonon propagator. The fermion field z a 
can be integrated out in 1/N expansion approximation, 
and the effective action of bosonic modes r, 8 can be ob- 
tained. The parameters r sa and A sa are determined by 
the saddle point condition, namely that the first order 
functional derivative of the effective action of bosonic 
modes vanishes. The Gaussian fluctuations with re- 
spect to this saddle point configuration are described by 
Sr = r — r sa and «r sa O. These Gaussian fluctuations play 
important role for impurity susceptibilities. 18 

In analogy with the Kondo spin singlet ground state 
we look for the orbital pseudospin singlet state, which is 
characterized by the condition 

run • (9) 

The electron spectral function in saddle point approxi- 
mation is given by 

A V = l (e-i f y + A* > A = Aor - (10) 
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In saddle point approximation the coherent Kondo peak 
structure is captured, but the incoherent high energy fea- 
tures are missing. ^From Eq. I jlOjl we can identify A with 
the (renormalized) Kondo energy scale (thus r^ a -C 1) 
and e f with the position of Kondo resonance peak. A and 
if are determined by the following saddle point equations 
(atT = 0): 
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D is the energy cutoff of the order of bandwidth of con- 
duction electron, and C is a numerical constant of order 
one. Tl(iu) is the polarization function 



U(iu) = - 
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If the electron-phonon coupling g is put to zero, Eq. i|ll|) 
and Eq. i|12|) reduce to the mean-field equation for the 
Kondo singlet ground state obtained by Coleman, [see 
Eq. 2.34 of Ref. 0| The saddle point equations can be 
solved numerically, but here we are content with the ap- 
proximate solution for the special case q — 1/2 (N = 2). 
In this case Eq. (|12f> implies that |e/| <C A. This also im- 
plies ef ~ —A sa < 0. Now Eq. (jl 1|> can be approximated 
to 
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The effect of the interaction with JT phonon is encap- 
sulated in the second term of the above equation. Using 
the explicit form of H(iw) we obtain (?/ is a numerical 
constant of order one) 



A 



where 
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for £1 < A, 
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2 A^ 
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is the unrenormalized Kondo temperature of the infinite- 
U asymmetric Anderson model obtained by Haldane us- 
ing poor man's scaling! 2 ^ El is the lattice relaxation en- 
ergy defined by 



decreases the Kondo temperature, which agrees with the 
previous expectation. For AH model the Kondo temper- 
ature increases by the interaction with Holstein phonon. 
The Kondo temperature of the infinite-C/ AH model ob- 
tained with the slave boson method isi& 
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To refine our discussion let us specify the reasonable rel- 
ative energy scales: |e/| > A > ft > (A > El)- In this 
parameter regime the renormalization of Kondo tempera- 
ture in AJT model is much larger than that in AH model. 
Apart from the magnitude of renormalization the depen- 
dence on various parameters is also markedly different 
between two models, reflecting the essentially different 
renormalization mechanism. 

It is interesting to attempt to understand the decrease 
of Kondo temperature by JT interaction from alternative 
approach. Let us consider the antiadiabatic limit where 
the phonon frequency is the highest energy scale, which 
in our situation implies that > |e/| > Ao > (A > El). 
Clearly this limit does not have much physical relevance, 
but it does provide some mathematical simplifications. In 
the antiadiabatic limit, the bare phonon propagator be- 
comes local in time Dq(t—t') — > 8(t— t')/MQ, 2 . Now the 
integrated JT electron-phonon interaction also becomes 
local in time (note that /, y are Grassman numbers and 
are not operators, in Lagrangian formulation) 

-E L T z T z ^+2E L f\f ] flfi. 

The right-hand side of the above equation is of the same 
form and of the same sign as the local Coulomb repul- 
sion, so that it increases the local Coulomb repulsion. 
To address the renormalization by JT interaction in this 
setup we had better consider finite-U Anderson model. 
For simplicity let us take the symmetric Anderson model. 
^From the Kondo temperature of the symmetric Ander- 
son model 22 with renormalized U one can obtain 



1 
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(17) 



which is similar to our result Eq. I|14l) with the replace- 
ment of Ao — > fi. Even if the antiadiabatic limit is not 
physical, this result gives a supporting evidence for the 
validity of our result. 

The impurity susceptibilities can be computed by in- 
serting source fields to the effective action Eq. (J7J and by 
integrating out z a and all bosonic modest The orbital 
susceptibility is defined by 
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This energy scale is also often called the polaron en- 
ergy. We find that the JT electron-phonon interaction 



In the Gaussian approximation one can obtain (C is a 
numerical constant) 

XoM =N\- CTI(icj) + C 2 g 2 D {iu)Tl 2 (iLu)\ , (19) 
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where H(iui) is given in Eq. I|13|) . Even though the details 
are not presented here, it should be be noted that the 
fluctuations of bosonic modes (Sr, O) do not contribute 
to the orbital susceptibility. Evidently, Eq. (|19|l are the 
first two terms of random phase approxiamtion (RPA)- 
type expansion, so that for small g 2 we can write 



XoO'w) 



N 



-Cn(iw)]- 1 -g 2 D (iuj)' 



(20) 



At high frequency w > Q the renormalization by JT 
phonon becomes negligible. In the low frequency limit 

LO -> 



Xo{iu -> 0) 
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The result Eq. H21[) clearly demonstrates the JT inter- 
action strongly enhances the orbital susceptibility. This 
again reflects the fact that the JT phonons couples to the 
orbital pseudospin channel directly. In sharp contrast to 
this result, the spin susceptibility of AH model is not 
renormalized by the interaction with Holstein phonon to 
the same order of approximation^ 
The charge susceptibility is defined by 



Xc(t 



(Sn f (r)Sn f (r')), 



(22) 



where n f = J2a s l s a = T,a z l z a and Sn f = nj - (n f ). 
The fluctuations of bosonic modes play very important 
role for the charge susceptibility. In the Gaussian ap- 
proximation we find that in the limit to — > 0, for the AJT 
model, 



Xc(i^ 



(23) 



Apparently the JT interaction seems to have no effect on 
the charge susceptibility, but this is misleading since the 
renormalization by JT interaction is reflected in A [Eq. 
ifTIJl] . Thus Eq. signifies that the JT interaction sup- 
presses the charge susceptibility. On the other hand, the 
charge susceptibility at zero frequency of the AH model 
is given by 



Xc(i^ -> 0) 
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and we can see that the chare susceptibility is enhanced 
but by a very small amount^ 

In the framework of DMFT there exists a systematic 
procedure linking the local impurity susceptibility to the 
susceptibility of the corresponding lattice system*^ In 
view of the fact that the orbital excitations can be de- 
tected by inelastic x-ray scattering, our result on the or- 
bital susceptibility can have important implications for 
the inelastic x-ray scattering studies of orbital physics. 
The renormalization of Kondo temperature can be also 
experimentally checked, for example, by the isotope ef- 
fects experiments on the transition-metal oxide systems 
and multiferroics. 
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